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PREFACE 


This  report  presents  the  results  of  studies  conducted  by  the 
Department  of  Civil  Engineering,  Louisiana  State  University,  to  support 
the  development  of  instruments  for  measuring  waves  and  sediments  in  the 
coastal  zone.  The  research  was  conducted  in  cooperation  with  the  Naval 
Coastal  Systems  Center  and  with  Dynex  Consulting  Company. 

Funding  for  this  study  was  provided  by  NCSC  and  the  Office  of  Naval 
Research  under  ONR  Contract  No.  N00014-79-C-0342  and  their  support  is 
gratefully  acknowledged. 
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I.  INTRODUCTION 

This  report  presents  the  results  of  progress  in  the  development  of 
a  predictive  scheme  for  computing  breaker  characteristics  from  wave 
measurements  outside  the  surf  zone.  This  research  was  conducted  in 
conjunction  with  an  instrument  development  program  of  the  Naval  Coastal 
Systems  Center.  The  NCSC  program  is  aimed  at  improving  UDT/SEAL  recon- 
naisance  through  the  development  of  a  operational  instrument  package 
which  will  provide  insitu  data  needed  to  determine  the  "effective  surf" 
(1)  at  a  beach  site. 

The  objective  of  this  study  was  to  synthesize  existing  knowledge 
about  breakers  and  surf  and  develop  predictive  models  to  be  used  with 
the  NCSC  or  similar  instruments.  A  secondary  objective  was  to  determine 
the  usefulness  of  a  penetrometer  for  measuring  marine  sediment  prop¬ 
erties. 

The  NCSC  instrument  will  provide  a  time  history  of  hydrodynamic 
pressure  and  the  two  horizontal  components  of  the  current  at  a  point 
seaward  of  the  breaker  and  surf  zone.  This  type  of  data  represents  the 
boundary  condition  for  surf  zone  wave  processes.  The  problem  is  to 
develop  a  scheme  for  using  boundary  value  data  from  offshore  instruments 
with  the  aid  of  either  a  mini-computer  (HP  9825  or  HP  9830)  or  a  large 
computer  (IBM  3033)  to  compute  surf  zone  parameters  such  as: 

1)  Average  breaker  height  (H)  and  maximum  breaker  height  (H_.) . 

2)  Average  breaker  position  (X)  and  maximum  breaker  position 

<y 

3)  Breaker  type. 

4)  Average  longshore  current  (V^)  and  maximum  longshore  current 

(V* 


To  be  of  engineering  use,  these  parameters  should  be  predicted  within 
♦10%  for  breaking  waves  up  to  a  maximum  height  of  10  feet  and  +20%  for 
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II.  ANALYTIC  PREDICTION  MODEL 

The  wave  theory  used  for  the  analytic  predictive  model  is  linear 
wave  theory.  This  theory  is  able  to  treat  waves  in  intermediate  and 
shallow  water  and  is  easily  programmable  onto  a  mini-computer.  It  is. 
of  course,  limited  to  waves  of  small  height  to  wave  length  ratios  (low 
steepness)  and  small  wave  height  to  water  depth  ratios.  It  is  used  in 
this  predictive  scheme  as  an  initial  step  in  the  development  of  the 
model/instrument  system. 

It  is  assumed  that  a  general  bathymetric  chart  (or  a  UDT  nearshore 
hydrographic  survey  chart  (NHSC)  (2))  is  available  for  the  beach  site. 

A  second  assumption  is  that  the  instrument  is  placed  at  a  known  distance 
off  the  bottom  (d^)  and  that  the  current  meter  is  oriented  so  that  the 
component  axes  are  parallel  (y  direction)  and  perpendicular  (x  direc¬ 
tion)  to  the  shoreline  (see  Figure  1).  It  is  further  assumed  that  the 
instrument  outputs  are  sampled  at  a  constant  time  interval ,  At,  over  an 
interval  of  time,  T,  to  produce  N  data  samples,  where  N  equals  T/At. 
Usually  T  will  be  in  the  range  of  5  to  20  minutes  and  At  will  be  .23  to 
1  sec. 


1.  Preliminary  Computations 

The  preliminary  computations  include  determining  the  mean  pressure 
(P),  the  pressure  change  between  wave  crest  and  wave  trough  for  each 
wave  (AP ^ ,  j*l,,,,M),  the  wave  period  for  each  wave  (T^),  the  mean 
velocity  components  (V  and  U)  and  the  wave  crest  current  components  (U 


and  V  .) . 
cj 


The  wave  parameters  computed  for  the  location  of  the  instrument 


are: 


A.  BEACH  PROFILE 
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Figure  1 .  Schematic  diagram  of  beach  and  definition  of  variabl 
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water  depth 


h  = 


P_ 
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and  for  each  wave: 


wave  height  H,  -  —  ■  t#su(lTfc/t  J) 

wave  length  ^  S - ^  •  TAUltt  {lltlr.  /C  ,'  ) 

"'i  i-i 


wave  angle 


2.  Predictive  Equations 

The  prediction  of  the  surf  zone  parameters  involves  a  set  of 
equations  which  must  be  applied  simultaneously.  The  equations  use  the 
instrument  site  wave  parameters  IT  ,  T^,  L^,  and  to  predict  the  same 
parameters  at  the  break  point  «b,  Tb>  Lb,  ©b.  It  is  assumed  that  the 
period  of  each  wave  remains  constant  so  that  T..  *  Tfe ^ .  The  equations 
are  as  follows  (3): 

1)  Breaker  Depth  Kkj  3  (l  *  4.7  W  4  *2  H*)  M-*,S 


2)  Breaker  Angle 


,  Cbr-[<|K* 


Vi 


3)  Breaker  Wavelength  L  kj  ' 


4)  Breaker  Height  M  bj  -  *  Kfj  •  l^sj 
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The  set  of  breaker  values  are  computed  for  j  =  1,,,,M.  For  this 
the  maximum  breaker  height  can  be  identified,  furthermore 


=-  j 

-W-  j»i 


The  positions  of  1^  and 

The  simultaneous  solution  of  the  set  of  equations  needed  to  make  a 
prediction  is  too  complex  for  a  practical  scheme.  Rather,  the  solution 
can  be  approximated  by  an  interaction  scheme  in  which  initial  values  for 


hj^^  can  be  located  on  the  nearshore  chart. 


selected  parameters  are  computed  from  the  input  data.  The  formula  for 
the  initial  value  of  the  wave  height  is 


<  =  .*«  **r*  M  jWk 


With  this  initial  value,  the  set  of  equations  in  the  surf  predic¬ 
tion  model  can  be  computed  explicitly  in  order,  that  is,  finding  , 
use  this  to  find  an<*  finally  .  The  difference  between  and 
should  be  small  (a  few  percent)  and,  it  not,  can  be  used  in 
place  of  in  a  second  iteration. 


Breaker  Type  and  Longshore  Current 


The  prediction  of  breaker  type  is  based  upon  breaker  wave  charac¬ 
teristics  and  the  beach  slope.  The  empirical  formula  suggested  is: 

H  y  1,2.  spiiu** 

Lb* 

£  £  jj  Vi  L  I.Z  pLUMCjlUfr 

^  suRt,i«b  BeeAfcefc 

Ui-m 

The  prediction  of  longshore  currents  is  based  upon  the  theory  for 
linear  waves  on  a  plane  sloping  beach.  The  formula  for  the  maximum 
longshore  current  and  the  longshore  current  at  the  average  breaker 
point  are  given  by  (4) : 

Vy#  -  20.*7  ^vw  f^5b)  1  SIKi  2.6b 

and 

V  j>  -  *  (oO  Yw 

4.  Verification 

The  analytic  model  presented  in  this  report  provides  a  means  for 
making  predictions  of  the  desired  surf  parameters.  The  model  is  limited 
in  application  because  of  the  assumptions  contained  in  the  underlying 
wave  theory.  The  model  does  not  account  for  wave  induced  pressure  near 
the  break  point  changes  which  differ  from  linear  wave  theory,  resulting 
from  wave  crest  accelerations  (5).  Furthermore  the  bottom  topography  is 
assumed  to  be  planar  and  therefore  real  beach  profiles  and  longshore 
bars  cannot  be  treated.  Finally,  the  longshore  current  predictions 
neglect  temporal  variations  in  current  speed  and  spatial  variations 
associated  with  nearshore  topography. 


Verification  of  the  mode]  accuracy  has  been  attempted  using  the  NCSC 
instrument,  as  reported  at  the  Oceans  80  Conference  (see  Appendix  1  for 
reprint  of  article).  The  mode]  predictions  of  average  breaker  heights 
(Table  1  of  reprint)  were  about  .10  to  .15  meters  too  high  for  average 
breaker  height  of  about  1  m.  This  could  be  accounted  for  by  an  in¬ 
creased  frictioned  attenuation  coefficient  (Kp) . 

5.  References 
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III.  NUMERICAL  PREDICTION  MODEL 


1.  Introduction 

Theoretical  studies  of  gravity  waves  in  incompressible  fluids, 
begun  in  the  early  nineteenth  century,  have  yielded  predictive  formulae 
of  successively  increasing  accuracy.  All  of  these  analytical  triu^s 
suffer,  however,  from  the  common  limitation  of  referring  to  physically 
simple  geometric  situations,  in  particular  a  wavefront  moving  over  a 
plane  bottom  having  a  constant  depth  or  slope.  In  order  to  obtain  a 
solution  to  the  general  problem,  we  are  compelled  to  introduce  discre¬ 
tization  of  the  space  and  time  dimensions,  and  solve  the  governing 
equations  for  the  dependent  variables  at  all  necessary  points  of  the 
four-dimensional  space-time  grid.  This  means  converting  the  analytic 
equations  to  finite-difference  form. 

However,  the  numerical  problem  posed  by  the  general  finite- 
difference  equations  is  still  formidable  on  account  of  Its  spatially 
three-dimensional  character.  In  the  simplest  physical  case  of  a  homo¬ 
geneous,  isothermal,  and  incompressible  fluid,  there  are  four  dependent 
variables:  u,  v,  w,  p — representing,  respectively,  the  three  components 
of  velocity  parallel  to  the  coordinate  axes  x,  y,  z,  and  the  pressure — to 
be  calculated  at  every  internal  grid  point,  and  the  water  level  to  be 
calculated  at  points  in  the  horizontal  plane.  For  nonisothermal  or 
inhomogeneous  liquids,  two  more  dependent  variables,  the  temperature 
and  density  must  be  added.  In  all  applications  to  geographic  situ¬ 
ations  where  many  wavelengths  are  involved  (short  wave  propagation)  the 
computer  storage  required  is  likely  to  be  excessive.  (At  least  ten  grid 
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points  are  needed  to  define  a  wave  profile).  Thus  the  modeler  will 
probably  find  himself  obliged  to  shrink  the  boundaries  of  his  investiga¬ 
tion  to  an  insignificant  si2e. 

Given  a  viable  grid,  however,  there  still  remains  the  problem  of 
specif iying  realistic  open  water  boundary  values  for  all  dependent 
variables  as  a  function  of  time,  over  the  vertical  plane. 

A  possible  way  out  of  these  difficulties,  and  one  which  drastically 
simplifies  the  numerical  procedure,  is  to  integrate  the  equations  of 
motion  from  bottom  to  surface  and  so  remove  one  spatial  dimension  explicitly 
from  the  problem.  All  dependent  variables  then  need  only  be  solved  for 
in  the  x-y  plane. 

The  modeling  of  shallow  water;  or  long  waves  (conventionally  defined 
as  ones  for  which  the  wavelength  Exceeds  twenty  times  the  mean  depth)  by 
means  of  vertically  integrated  equations,  has  received  much  attention. 

Ever  since  Hansen  (1956)  introduced  a  viable  finite-difference  formulation 
of  these  equations,  numerous  authors  have  applied  his  technique  to 
various  seas,  lakes,  and  bays.  The  long  waves  of  interest  to  oceanographers 
are  principally  those  caused  by  the  sun  and  moon  (tides),  and  those 
arising  from  strong  winds  and  atmospheric  fluctuations  (storm  surges). 
Successful  predictions  of  these  phenomena  have  been  obtained  using  the 
Hansen  equations  solved  explicitly,  implicitly,  or  a  mixture  of  both. 

A  desirable  first  approach  to  the  solution  of  the  general  wave 
problem  in  real  geographic  situations  would  be  to  apply  the  numerical 
techniques  already  proven  for  solving  the  long  wave  equations.  The 
"general”  wave  equations  would,  however,  need  to  retain  a  higher  order 
of  validity.  In  obtaining  soluble  long-wave  equations  it  has  been 
customary  to  make  two  simplifying  assumptions:  (1)  that  the  horizontal 


'  / 
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velocity  vector  Is  constant  from  bottom  to  surface,  and  (2)  that  the 
vertical  acceleration  is  so  small  that  the  hydrostatic  approximation  is 
a  good  one.  In  the  general  case  these  restrictions  would  have  to  be 
relaxed . 

To  account  for  the  unknown  variability  of  the  horizontal  velocity 
component  in  the  vertical  direction,  and  a  finite  vertical  acceleration, 
while  keeping  the  problem  amenable  to  a  Hansen-type  approach,  requires 
the  help  of  analytical,  or  closed-form  solutions.  These  necessarily 
refer  to  the  simplest  geometric  situation;  but  used  in  conjunction  with 
the  vertically  integrated  Navier-Stokes  equations  it  is  expected  that 
they  will  allow  a  good  first  approximation  to  the  solution  of  the  general 
problem. 

* 

The  recognition  of  a  non-zero  vertical  acceleration  in  the  fluid 
flow  is  tantamount  to  accepting  a  pressure  distribution  in  the  vertical 
deviating  from  the  hydrostatic.  Conversely,  with  a  pressure  distribution 
in  the  vertical  taken  from  analytic  or  "zero  order"  theory,  the  vertical 
acceleration  is  implicitly  accounted  for. 

2.  Derivation  of  the  Vertically  Integrated  Equations  for  the  General 
Wave  Problem 

The  three-dimensional,  ensemble-averaged  equations  of  motion  for  a 
homogeneous,  incompressible  fluid,  with  turbulent  stresses  on  vertical 
planes  and  molecular  viscous  stresses  neglected,  can  be  written 


—  +  9v  +  ^w  *  0 
3x  3y  Tz 


3u  +  3u2  +  3uv  +  3uw  -  2 ftzv 

3t  3x  3y  3z 

3v  3vu  +  3v?  +  3w  +  2ftzu 

3t  3x  3y  3z 


"  .li.  l£  "  i  ^xz 

0  3x  p  3z 

■  -1  3p  -  1  tVtyz 

D  3y  p  3z 


3w  3wu  3wv 

3c  3x  3y 


+  2(f^v  -  J^u) 


1  3p  1  3^zz 

p  ®  '  p 


(1) 
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where  fl  is  the  earth's  angular  velocity  vector,  rxz»  TyZ*  Tzz  *re 
turbulent  stress  components  on  the  horizontal  plane  in  the  directions  z, 
y,  z,  respectively,  g  is  the  acceleration  due  to  gravity;  and  other 
symbols  have  been  defined.  The  Coriolis  components  of  acceleration  have 
been  simplified  by  omission  of  terms  Involving  w,  on  the  assumption  that 
w  «  u,  v. 

We  eliminate  the  third  momentum  equation  at  the  outset  by  adopting 
the  pressure  distribution  from  small  amplitude,  frictionless  wave  theory 

(see,  for  example  Zppen  (1966)): 

p  -  -gpz  +  gpcosh[k(h  +  z)l.C  o) 

cosh[kh] 

In  this  equation  the  external,  or  atmospheric  pressure  is  implic¬ 
itly  zero,  and  the  origin  of  z  (measured  positively  upward)  is  located 

* 

at  the  still  water  surface,  h  is  the  depth  of  the  bottom  below  this 
surface.  The  instantaneous  height  of  the  liquid  surface  relative  to  the 
plane  z  *  0,  is  Z.  k  is  Che  wave  number,  defined  by 

k  -  2*/L 

where  L  is  the  wavelength.  For  the  purpose  of  the  ensuing  analysis  both 
k  and  h  are  considered  -  in  contradistinction  to  linear  theory  -  to  be 
functions  of  x  and  y.  Figure  1  illustrates  the  definitions  of  z,  h  and  s. 


(  See  Figures  Section  of  Report  ) 

Fig.  1  Definition  Diagram  of  z,  h,  and  £. 

In  the  derivation  of  equation  (2)  the  surface  boundary  condition 
was  applied  not  at  the  actual  surface  z  -  C,  but  at  z  »  0.  Conse¬ 
quently,  the  equation  has  validity  only  from  z  *  -h  to  z  «  0.  For 
positive  z  one  must  assume  a  hydrostatic  pressure  distribution. 
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Write  iff(z)  ■  eosh[k(h  4-  z)]  -h  £  z  £  0  (3) 

cosh[kh] 


*P  has  been  termed  the  "pressure  response  factor".  From  equations  (2) 
and  (3) 


13£«g<J;3C,13o  =  g^3C 

p  3x  3x  p  3y  3y 


With  the  last  two  relations,  and  replacing  2J1  by  the  Coriolis 

z 

parameter  f,  the  equations  to  be  solved  become 

3u  +  3v  +  3w 
3x  3y  3z 

3u  +  3u2  +  3uy  +■  3uw  -  f v  *  -gt[>  _3^_  -  1^  3*^cz 

3t  3x  3y  3z  3x  p  3z 

3v  +  3vu  +  3v2  +  3vw  +  fu  ■  -gty  3^  -  1^  3~^yz 

3t  3x  3y  3z'  3y  P  3z 


Equations  (4)  must  be  integrated  over  z  between  the  limits  -  h  and 
£.  Consider  firstly  the  integration  of  Here,  in  accordance  with  the 
previous  remarks,  the  upper  limit  must  be  changed  to  z  -  0.  The  error 
will  not  be  serious  if  the  small  amplitude  condition:  |c|«  h  ,  holds. 
From  equation  (3)  0 

/  ip  dz  ■  1  tanh[kh]  (s% 


Now,  another  result  of  the  linear,  or  small  amplitude  theory  is 


that 


a2  -  gk  tanhfkh] 


(6) 


where  0  is  the  angular  frequency  of  the  wave,  a  constant.  Eliminating 
tanh  [kh]  from  equation  (5)  using  equation  (6),  we  then  obtain  approxi¬ 
mately  r  _ 

r  gip  3£  dz  ■  a2  3c 

-h  lx  k*  lx 


1  / 


(7) 
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By  definition  of  the  phase  speed,  S 

a2  -  s2  <8> 

F 

For  long  waves,  when  <l>  ■  1  is  a  good  approximation,  the  right  hand 
side  of  equation  (7)  is  replaced  by 
g(C  +  h) 

3x 

Boundary  conditions  that  must  be  employed  in  the  vertical 
Integrations  are: 

3C  +  u3c  +  v3g  •  w  z  m  C 

3t  3x  3y  (9) 


3h  +  u3h  +  v  3h  *  -w  z  *  -h  (10) 

3t  3x  3y 


Using  Leibnitz's  Rule  with  conditions  (9)  and  (10)  and  relations  of 
type  (7),  we  vertically  integrate  equations  (4)  between  the  limits  z  »  -h 

and  z  -  5  to  yield 

3£  +  3_  /C  udz  +  3_  vdz  »  0 
3t  3x  -h  3y  -h 

3_  udz  +  3  u2dz  +  3_  uvdz  -  f/^vdz  *  -<£  +  t*  -  Tb 
3t  -h  cix  -h  3y  -h  -h  -I?  3x  —  * 

3  vdz  +  vudz  +  v2dz  +  f/^  udz  *  -o2  jte  +  xs  - 

3?  -h  'h  3y  ~h  -h  F  3y  -2 - X 

P 


s  b 

and  T  are  the  x-components  of  the  surface  and  bottom  stresses, 
respectively.  Similarly  with  the  y-components. 


We  now  follow  Leendertse  (1967)  in  writing 
u  -  U(1  +  eu) 
v  -  7(1  + 


i 


(12) 


where 


V 
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U  -  1  z5  udz  , 

h  +  C  -h 


Z5 


vdz 


h  +  Z  -h 

Thus  U  (x,  y)  aod  V  (x,  y)  are  the  vertically  averaged  values  of  u 
(*»  7.  *)  and  v  (x,  y,  z),  respectively;  and  (x,  y,  z)  and  Ev(x,  y, 
z)  express  the  fractional  deviations  of  u  and  v  from  U  and  Vt  respec¬ 
tively,  at  x,  y,  z. 

Define  the  total  height  of  the  water  column 

H  -  h  ♦  C 


(13) 


Then,  using  equations  (12)  and  (13),  the  left  hand  sides  of  equa¬ 
tions  (11)  become 

3£  +  3HU-+  3HV 
3t  3x  Zy 

3HU  +  3_  {U2(H  +  Z?e2  dz)}  +  {UV(H  +  ZCe  e  dz)}  -  fHV 
3t  3x  -h  u  3y  -h 

3HV  +  3_  {VU(H  +  J’ceuevdz)}+  3_  {V2(H  +  Z^ejdz)  }  +  fHU 
3t  3x  3y  -h 


Assume  isotropy  of  the  horizontal  velocity  stratification.  That  is 

e  (*»  7.  *)  ■  e  (*,  y»  *)  *  e(*.  y»  *)  (14) 

U  V 

and  let 

H  +  Z?  e2dz  -  aH  (15) 

-h 


a(x,  y*  O  may  be  termed  the  "local  velocity  stratification  factor" 
at  time  t.  Always 

a  >  l 

With  assumption  (14)  and  definition  (15),  and  after  some  reduction. 


i 


_ _  i 
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Che  equations  governing  the  flow  become 

3£  +  3HU  +  3HV  -  0 
7t  3x  3y 


3U  +  aD  3U  +  aV3U 
3t  3x  Ty  " 


fV 


-az  3;  +  t8 

1?H^  x 


-  t: 


oh 


(16) 

{D23a  +  UV  3a  +  0  (1  -  a  3c} 

H  3?  I  Jt 


3V  +  aD  3V  +  aV3v  +  fU 
3t  3x  ijy 


-a2  3c 


{VU  3a  +  V2  3a  +  V  (1  -  a)3c> 
3^  Hy  H  3F 


2  2 

If  in  equations  (16)  we  put  a  -  1  and  replace  o  /k  H  by  g,  we 

obtain  the  Ha  nsen  long  wave  equations.  From  (6)  we  see  that  the  second 

replacement  amounts  to  the  equality 

tanh  [kh]  -  kR 
s 

This  will  be  closely  approached  for  shallow  water  (or  long)  waves 
provided  Id  <K  h. 

3.  Evaluation  of  the  Velocity  Stratification  Factor 

It  is  desired  to  evaluate  a  in  closed  fern,  using  the  results  of 
linear  wave  theory. 

On  the  same  level  of  approximation  as  equation  (2) 

u  *  jjk  cosh[k(fi  +  z)1 .  C  -h  ^  z  £  0  (17) 

<3  cosh[kh] 

From  equations  (3),  (5)  and  (17)  we  have 

f^udz  =*  gCtanhfkh] 

-h  a 

Hence,  approximately,  by  definition  (13)  and  the  last  result 

U  ■  gCtank[kh]  (18) 

<JH 

and  therefore,  approximately 

C  -  UaH 

gtanh[kh] 

Substituting  the  last  approximation  for  5  in  (17)  yields 

u  ■  UkH  cosh[k(h  +  z)1 
sinh[kh] 
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whence,  comparing  with  the  first  of  equations  (12) 

1  +  e  -  kH  cosh[k(h  +  z)  1  (19) 

siohfkh] 

Now,  by  definition  of  £,  we  have 

/C  (1  +  e)2  dz  =  H  +  e2  dz 
-h  -h 

Hence,  from  equations  (IS),  (19),  and  the  last  result 

ctH  -  /?  k2H2  cosh2fk(h  +  z)]  dz 
-h  sinh2[kh] 

from  which  2 

a  -  1  (kH)  (1  +  sinh(2kH]) 

2  (sinhlkh])2  2  kH 

It  v*il  be  desirable,  for  computational  reasons,  that  a  should  be  a 
function  of  position  only.  To  remove  the  time  dependency  of  a ,  we 

substitute  h  for  H  in  the  last  result,  and  so  obtain 

3  "  JE  (kh)2  ,  (1  +  sinhKkh]) 

2  (sinh[khj)  2  kh  U0' 

* 

Equation  (20)  again  only  holds  with  the  small  amplitude  assumption. 
Some  computed  values  of  a  between  the  asymptotic  limits  1  and  kh/2 
are  shown  in  Table  1. 

Table  1 

Some  Computed  Values  of  a 
(  See  Tables  Section  of  Report  ) 


It  appears  that  the  Hansen  equations  are  adequate  well  into  the 
"intermediate  wave"  zone,  and  should  yield  fair  accuracy  for  wavelengths 
as  short  as  3h. 

4.  The  Two-Dimensional  Numerical  Model 

Prior  to  the  present  study  a  two-dimensional,  vertically -averaged, 
long  wave  computer  model  had  been  developed  by  McHugh  (1976),  utilizing 


i 


a  mixed  implicit-explicit  method  of  solving  the  finite-difference 
equations  derived  from  Leeadertse  (1967).  The  discretization  was  that 
shown  in  Figure  2. 


Fig.  2  Discretization  of  C,  h,  U,  V 
A  full  account  of  the  genesis  of  the  model,  which  was  designed  to 
accomodate  stretching  of  the  grid  lines  and  a  time-varying  land/water 
boundary,  is  given  in  McHugh  (1976).  The  model  was  subsequently  modi¬ 
fied  to  allow  randomly  spaced  gri^  lines,  and  this  is  the  form  utilized 
in  the  present  work. 

In  order  to  generalize  the  long  wave  model  to  simulate  waves  of 
arbitrary  length,  it  was  merely  necessary  to  insert  into  the  momentum 
equations  the  factors  ®  and  the  correction  terms  in  braces  as  shown  in 

equations  (16);  also  to  replace  g  in  the  pressure  driving  terms  by 

2  2  2 
cr/k  H  -  S*/H. 

■Hie  bottom  stress  was  represented  through  the  Ch&zy  coefficient, 

C(x,y,t) .  Thus,  T*  »  pgO(U2  +  V2)1/2  ,  Tb  »  pgV(u2  +  V2)1/2 

C2  7  C2 

C  values  were  computed  from  Manning's  formula: 

C  -  1.486  HI/6  (H  in  ft.) 

n 

The  resultant  finite  difference  equations,  which  are  very  involved, 
have  been  presented  as  Appendix  1.  Referring  to  this  appendix,  the 
superscript  n  (1,  2,  3...)  denotes  the  time  level  n&t.  Barred  quan¬ 
tities  with  subscript  u  or  v  denote  two-point  interpolated  values  at  the 


positions  of  0  or  V  indicated.  Quantities  V  *  .  and  0  *  denote  four- 

*  *j  *  »j 

point  interpolations  of  V  and  0,  respectively,  centered  at  the  locations 
of  j  and  y  All  of  these  interpolations  take  into  account  the 
possible  non-uniform  grid  spacing. 

The  quantities  1^,  etc.  are  multipliers  for  determining  deriv¬ 
atives  when  the  grid  spacing  is  irregular  or  stretched.  In  equation 
A(i)  for  instance,  for  irregular  spacing: 

T  I  1 


T  ' 

i+1/2 


*i+l  -  xt 


The  notation  derives  from  the  form  when  I  ■  I(x)  and  J  ■  J(y)  are 

stretching  functions  such  that 

3__3_dI.jL.3_dJ 
3x  91  dx  3y  3J  dy 

or,  in  finite  difference  operator  form,  writing  I*  for  dl/dx  and  J'  for 


dJ/dy 


A__«  A_  I*  ,  A  ,  A_  J’ 
Ax  AI  Ey  m  AJ 


Lover  case  subscripts  have  been  used  in  the  equations  for 
convenience . 

Quantities  S.  .  and  O  were  calculated  at  the  centers  of  the 
L»J  I*J 

corresponding  grid  squares.  Two  sorts  of  derivatives  of  a  were  used: 
those  centered  on  the  £ -points,  and  those  centered  on  the  D  and  V-points. 
The  latter  are  distinguished  by  additional  subscripts  u  or  v  in  the 
manner  already  discussed.  Thus  for  example 

°*ifJ  *  <>1+1,3  ~  °i-l» J 

l/2(Xi+^  +  xi  xi-l 

»xuitj  *  «i,J  ~  <*1-1, j 

l/2(x1+i  -  xi-i)  _ 

with  analogous  expressions  for  Oy.  and  ayv.  . 

For  the  constant  boundary  model,  k(x,  y) ,  S  (x,  y) ,  and  a(x,  y) 

are  calculated  once  and  for  all  time  over  the  entire  wet  field  from 

equations  (6),  (8)  and  (20),  respectively.  Newton's  iterative  method 


1 


1 


20 


vas  found  adequate  for  extracting  k  from  equation  (6).  The  first  guess 

of  k  was  the  deep  water  value,  kQ,  given  as  input  to  the  computer  program. 

When  the  wet  boundary  changes  with  time  (as  for  instance  on  a  beach) 

there  is  the  necessity  of  defining  S  and  a  in  the  newly  flooded  region. 

The  strategy  in  the  present  study  was  to  assume  a  default  value  a  ■  1, 

2 

and  a  default  value  S  /H  *>  g  in  the  initially  dry  regions  (h  ^0). 

The  implicit-explicit,  double  sweep  solution  technique  has  been 
described  by  Leendertse  (1967)  and  McHugh  (1976).  Its  principal  stages 
are:  (I)  Solve  implicitly  for  +  *  and  11°  +  *  along  each  grid  row. 

(ii)  Solve  explicitly  for  Vn  +  *  in  each  square  of  the  wet  field,  (iii) 
Solve  implicitly  for  CQ  +  ^  and  Vn  +  ^  along  each  grid  column,  (iv) 

Solve  explicitly  for  U°  +  ^  in  each  square  of  the  wet  field.  Repeat  the 
procedure  for  the  next  pair  of  time  steps,  and  so  on.  If  the  land/watcr 
boundary  varies  in  time,  this  must  be  established  periodically — eg., 
every  two  time  steps — and  allowance  made  for  volume  transport  into 
previously  dry  areas  (see  Appendix  2) . 

A  principle  of  the  computation  is  to  ignore  derivatives  when  they 
cannot  be  computed  from  the  central-difference  formula  owing  to  the 
proximity  of  a  land  boundary.  This  is  found  (Leendertse  (1967))  to  be 
beneficial  for  numerical  stability. 

Finally,  the  necessary  boundary  conditions  required  on  a  row  or 

column  are  either  the  water  level  C,  or  the  input/output  velocity  0  or  V 

at  all  times.  At  a  land  boundary  0  ■  V  »  0. 

At  the  start  of  the  model  it  is  customary  to  assume  C  »  O^cvri  U- V-  0 


everywhere  inside  the  wet  field. 

5.  Preliminary  Tests  of  the  Model 
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A  5  x  121  grid  was  selected  initially  with  step  lengths  Ax  ■  Ay  *  16  ft 
(4.38  m) .  Three  sides  were  closed  thus  forming  a  channel  with  dimensions 
64  x  1920  ft.  (19.52  x  585.60  m) .  For  this  geometry  it  will  be  convenient 
to  consider  the  x-axis  as  running  parallel  to  the  longer  dimension. 

Thus  U  is  the  onshore  velocity  component.  For  all  tests  with  the  channel 
a  standard  sinusoidal  disturbance  of  water  level  was  input  at  the  open 
boundary  having  period  8  seconds  and  amplitude  1  ft.  (0.305  m) .  The 
deep  water  wavelength  L^,  deduced  from  equation  (6)  with  tanh  [kh]  *  1, 
is  then  approximately  320  ft  (97.60  m) .  With  this  Lq  the  16  ft.  grid 
line  separation  yields  20  grid  intervals  per  wavelength  in  deep  water. 

If  the  bottom  rises  the  wavelength  will  diminish,  and  it  may  be  desirable 
to  increase  the  grid  resolution  along  x  in  the  shallow  region.  Various 
resolution  schemes  were  considered,  maintaining  however,  a  uniform  grid 
interval  along  y  of  16  ft. 

Four  depth  contours  were  tested:  (a)  h  constant  at  150  ft.  (45.75 
m)  ,  (b)  h  varying  linearly  from  150  to  10  ft.  (45.75  -  3.05  m) ,  (c) 
varying  linearly  from  150  to  -  3  ft.  (45.75- -0.92  m) ,  and  (d)  h  varying 
linearly  from  150  to  -  3  ft.  except  in  the  region  of  a  submerged  bar. 
Contour  (a)  was  used  with  both  a  closed  (vertical  wall)  and  open  boundary 
"downstream”.  In  (b)  waves  reflect  off  a  vertical  wall;  in  (c)  and 
(d)  waves  reflect  off  a  beach.  Case  (b),  which  is  intermediate  in 
effects  between  (a)  and  (c) ,  may  be  called  the  "quasi-beach"  conf ig-  - 
uration;  it  was  chosen  to  eliminate  possible  numerical  distrubances 
arising  from  the  changing  land/water  boundary,  while  preserving  the 
effect  of  shoaling. 
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As  earlier  seated,  in  any  wave  computation  equations  (6)  and  (20) 
must  be  solved  at  all  positive  values  of  h  for  k  and  a.  Using  equation 
(8)  and  the  definition  L  ■  2ir/k,  the  wavelength  L  and  phase  speed  S  may 
be  obtained  for  every  grid  square.  Table  2  lists  some  values  of  L,  S, 
and  a  for  h-values  within  the  range  of  the  tests. 

Table  2 

Some  values  of  L,  S  and  a  as  Functions  of  h 


A.  Constant  h?  Deep  Water  Configuration 

With  h  -  150  ft.  it  may  be  deduced  from  Table  2  that  the  wave 
should  traverse  the  length  of  the  channel  in  approximately  47  seconds. 

The  computation  was  therefore  allowed  to  simulate  52  seconds,  with 
tabular  output  of  c,  U,  and  V  every  0.5  seconds.  A  time  step  was  initially 
selected  of  0.05  seconds.  This  is  well  below  the  Friedrichs-Lewy-Courant 
stability  limit  of  0.167  seconds  for  the  considered  grid  interval  and 
depth.  Results  for  £  and  U  with  At  *  0.025  seconds  agreed  in  two  places 
of  decimals,  hence  we  are  close  to  solution  convergence.  The  appropriate 
roughness  coefficient  for  short  waves  over  a  150  ft. -deep  mud  bottom 
being  unknown,  the  value  n  *  0.026  was  arbitrarily  chosen.  This  value 
is  typical  of  tidal  models. 


Fig.  3  Wave  Profile  for  Deep  Water  Configuration  (h-150ft.) 
£X-16ft>,  &t»0.05sec.,  n-0.026 


In  Figure  3,  which  ehove  the  wave  profile  at  48  seconds  (&(i,  2), 
1-1,.. 120)  It  can  be  seen  that  the  solution  is  numerically  well  behaved. 
The  solution  at  72  seconds  (not  shown)  was  also  smooth  and  noise-free. 

Figure  4  shows  U  plotted  as  a  function  of  t  for  the  grid  point  (80, 
2).  This  point  was  chosen  for  examination  as  being  In  a  region  where 
the  equilibrium  wave  profile  was  beginning  to  be  established,  yet  not  so 
close  to  the  solid  boundary  as  to  Invite  Interference  from  reflection 
effects. 


Fig.  4  U  as  a  Function  of  t  for  the  Deep  Water 
Configuration  at  i-polnt  80 
Note:  lft/sec  -  0.305m/sec. 

According  to  equation  (18)  ve  may  write  appr> wlmately  for  deep 

water 

U(z,t)  - 

oh 

Since  linear  theory  requires  £  to  have  the  form 

£  -  a  sin  (lex  -  <rt) 

we  have  for  the  amplitude  of  U,  with  a  -I  Ft  and  d*  -  2^/8  r acJ/stt 

O'  -  ag  -  0.273  ft/sec. 
oh 

The  maxima  and  minima  of  Figure 4  are  0.283,  -0.270,  0.294,  -0.262, 
0.290  ft  /sec  with  mean  absolute  value  0.280  ft  /sec.  The  difference 
from  O'  calculated  above  Is  less  than  2Z.  Since  the  maxima  and  minima 
appear  to  be  decreasing  In  absolute  value  the  agreement  would  probably 
have  been  even  better  for  a  longer  channel. 
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A  second  quantitative  comparison  between  the  model  and  theory  is 
afforded  by  the  so-called  Stokes  drift  velocity,  u*,  integrated  over  the 
depth.  Kinsman  (1965)  derives  u*  from  the  velocity  potential  in  the  form 

u*  -  a2k2  Se2kz 

whence 


u*  «  a2k2S  e2kz  dz  -  a2kS 
h  -h  2h 


With  S  ■  o/k  we  have 

tt*  *  a2c  *  0.00262  ft/sec. 

2h 

From  the  tabulated  output  of  U  at  i-point  80  we  estimate  numerically 
0  -  /TU(h  +  Odt  -  r  n  .  H.  At’ 

o  hT  3  'i 

h  Zj  AtT 

where  T  is  the  period,  Uj ,  are  tabulated  values  of  U  and  H  at  time 
t j ,  end  At*  is  the  constant  time  interval  (0.5  seconds).  The  j-summation 
is  taken  over  a  number  of  complete  cycles  of  0.  With  the  time  limits 
for  integration  t  ■  28.5  seconds  (U  •  0.000  ft. /sec)  and  t  ■  47.0 
seconds  (U  «  -0.010  ft  /sec  ) ,  we  have  ■  14.251  ft.  /sec.,  E^t* 

*  18.5  seconds,  and 


"0  -  0.00257  ft. /sec 

The  difference  between  U  and  u*  is  less  than  21, 

B.  Linearly  Varying  h:  Quasi-Beach  Configuration 

Figure  5  shows  wave  profiles  at  60  and  61  seconds  for  a  high 
resolution  grid  having  Ax  ■  16  ft  for  the  first  90  intervals;  &x  •  8 
ft.  for  the  next  20  intervals;  Ax  -  4  ft  for  the  next  40  intervals;  and 
Ax  -  2  ft  (0.61  m)  for  the  last  80  intervals.  The  total  length  of  the 


grid  is  the  same  as  in  subsection  A. 


Fig.  5  Quasi-Beach  Have  Profiles  at  60  and  61  Seconds. 

Ax»16,  8,  4,  2ft  ,  A  t**0.0125sec  . ,  n-0.026 

One  sees  clearly  a  breaking  wave  approaching  the  end  wall. 

In  Figures  6a  and  6b  the  time  has  been  increased  to  72  seconds.  The 

higher  resolution  grid  has  Ax  “  16  ft  for  the  first  90  intervals; 

Ax  *  8  ft,  for  the  next  20  intervals;  and  Ax  -  4  ft.  for  the  last  80 

intervals,  keeping  the  total  length  the  same  as  in  subsection  A. 


Fig.  6  Quasi-Beach  Wave  Profiles  at  72  seconds;  At  *  0.0125  sec  ,  n 
•  0.026  (a)  Ax  -  16  ft  ,  (b)  Ax  -  16,  8,  4  ft. 

Reflection  of  the  breaking  wave  at  the  wall  has  apparently  gener¬ 
ated  numerical  noise  propagating  back  into  the  channel.  This  noise 
exhibits  wavelengths  equal  to  simple  multiples  of  the  local  grid  spac¬ 
ing.  As  grid  resolution  increases, more  spurious  waves  are  generated  per 
unit  length  of  the  channel  in  the  noise  affected  area;  and  certain 
amplitudes  may  increase  so  much  as  to  cause  numerical  Instability,  which 
is  evidenced  by  fluctuating  non-zero  V-components  in  the  transverse 
direction,  and  transverse  oscillations  of  £. 


Elimination  of  Grid-Induced  Noise  by  Artificial  Smoothing 


An  approach  to  the  elimination  of  noise  which  has  been  found  useful 

by  tidal  modelers  is  the  so-called  "artificial  viscosity"  strategy.  One 

replaces  certain  dependent  variables  X°  .  in  the  finite  difference 

i»j 

equations  evaluated  at  the  old  time  level  n,  with  quantities  ]  defined, 

for  a  uniform  grid,  by 


• «!, j +  i/4<i  -  s)  (*;+lij  + 


+  x 


+•  xn. 


i-u-  i.J+1  i.J-1 


,) 


0  is  a  dimensionless  number  with  value  usually  between  0.99  and 
1.0.  A  value  of  0  less  than  unity  has  the  effect  of  lowering  peaks  and 
raising  troughs  relative  to  the  mean  of  the  variable  X°  hence  smoothing 
it.  In  the  finite  difference  equations  of  Appendix  1,  artificially 
smoothed  variables  are  enclosed  by  brackets  [  ]. 

Some  smoothing,  and  therefore  de-energizing  of  the  physically  real 
wave  will  also  result  from  this  ploy;  but  it  must  be  borne  in  mind  that 
real  world  waves  possess  internal  friction  which  the  model  represented 
by  equations  (16)  does  not.  Hence  one  may  be  able  to  live  with  a  slight 
smoothing  of  the  longer  wavelengths. 

The  16-8-4  grid  scheme  was  used  with  low  bottom  friction  (n  • 

0.026)  to  test  the  effect  of  introducing  0  <  1  into  the  finite-difference 
equations.  Figure  7  shows  wave  profiles  at  72  and  73  seconds  for  a 
number  of  values  of  6.  The  uppermost  trace  corresponds  to  Figure  6b. 


Fig.  7  Quasi— Beach  Wave  Profiles  at  72  and  73  Seconds  With  Artificial 
Smoothing.  Ax»l6,  8,  Aft  ,  At*0.0l25sec  ,  n*0.026 


Almost  Identical  traces  to  those  for  B  *  0.995,  At  *  0.0125  sec  , 
were  obtained  with  6  -  0.990,  and  At  *  0.025  sec  ,  showing  that  for 
B  <  1  the  solutions  become  more  sensitive  to  the  time  step  than  with 
6*1.  One  thus  has  two  degrees  of  parametric  freedom,  B  and  At,  to 
achieve  the  same  degree  of  smoothing. 

In  the  last  pair  of  traces  in  Figure  7  the  wave  breaking  effect 
appears  to  have  been  virtually  eliminated.  This  may  be  useful  for  some 
applications . 

Another  way  of  reducing  the  energy  of  numerical  turbulence  is  to 
introduce  substantial  bottom  friction  in  the  zone  of  wave  breaking. 

Figure  8  shows  wave  profiles  at  72  and  73  seconds  for  the  highest 
resolution  grid  (16-8-4-2)  and  a  Manning's  n  distributed  as  in  the 
figure  caption.  The  width  of  the  zone  with  n  *  5  is  approximately  54 
ft  (16.5  m).  The  unsmoothed  profiles  (B  *  1)  have  been  given  for 
comparison  with  smoothed  profiles,  B  *  0.997. 


Fig.  8  Quasi-Beach  Wave  Profiles  at  72  and  73  Seconds  With  Combined 
High  Bottom  Friction  and  Artificial  Smoothing. 

Ax  *  16,8,4,2  ft  ,  At  *  0.025  sec  ,  n  *  0.026,  h  ^  20  ft. ;  n  * 

0.2,  14  ft  <  h  <  20  ft  ,*  n  -  5,  h  <  14  ft. 

C.  Linearly  Varying  hi  Beach  Configuration 

This  configuration  involves  a  time-varying  land/water  boundary. 

For  a  brief  description  of  the  handling  of  this  problem  see  Appendix  2. 

Figure  9  shows  the  wave  profile  at  76  seconds  for  the  uniform 

coarse  grid  and  the  lowest  bottom  friction  without  numerical  smoothing. 


Grid-gene rated  noise  is  present  near  the  shoreline.  Increasing  the  28 
resolution  of  the  grid  for  this  low  bottom  friction  resulted  in  unstable 
solutions  for  all  grid  schemes  considered. 


Fig.  9  Beach  Vave  Profile  at  76  Seconds 

Ax  •  16  ft  ,  it  ■  0.0125  sec  ,  n  ■  0.026 

Using  the  16-8-4  grid  scheme  and  n  «  0.026,  all  values  of  g  down  to 
0.99  did  not  succeed  in  removing  the  transverse  numerical  instability  at 
76  seconds.  At  8  “  0.95,  the  wave  height  in  the  zone  x  >  1440  ft.  was 
less  than  0.2  ft  (0.06  m) ;  hence*  n  -  0.026  is  not  a  viable  Manning 
coefficient  for  the  beach  configuration  even  with  artificial  energy 
loss. 

In  Figure  10  we  have  wave  profiles  at  76  seconds  for  the  16-8-4 
grid  scheme  and  a  Manning's  n  distribution  as  given  in  the  caption.  The 
slight  noise  in  the  unsmoothed  trace  has  been  almost  eliminated  in  the 
trace  with  8  »  0.997. 


Fig.  10  Beach  Wave  Profiles  at  76  Seconds 

Ax  -  16,8,4,  At  -  0.025  sec  ,  n  «  0.026,  h  >  20  ft  ;  n  ■  0.2,  10 
ft  <  h  <  20  ft;  n  0.5,  h  <  10  ft.  “ 
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D.  Linearly  Varying  h-— Beach  Plus  Offshore  Bar 

Figure  11  shows  the  wave  profile  at  60  seconds  that  was  obtained 
for  the  bar  geometry  shown  inset,  using  the  uniform  coarse  grid.  The 
bar  ridge,  48  ft  (14.6  m)  wide  is  approximately  6  ft  (1.83  m)  deep  and 
368  ft  (112.2  m)  offshore. 


Fig.  11  B each-Wi th-Bar  Wave  Profile  at  60  Seconds 
Ax  ■  16  ft  ,  At  0.05  sec  ,  n  »  0.026 

There  is  probably  noise  present  in  the  wave  trace,  as  suggested  by 
the  experience  of  Figure  9.  Significant  results,  however,  could  doubt¬ 
less  be  obtained  by  a  judicious  choice  of  bottom  friction  over  the  bar, 
a  better  grid  resolution  in  the  same  neighborhood,  and  some  artificial 
smoothing. 

E.  Constant  hi  Deep  Water  Configuration  With  Open  Boundary  Downstream 
If  a  beach  or  other  kind  of  closed  downstream  boundary  is  not 
desired,  an  open  boundary  or  reactionless  gate  to  the  passage  of  the 
waves  may  be  simulated  by  an  appropriate  adjustment  of  either  the  water 
level  in  the  last  grid  square,  or  the  0-value  on  the  boundary,  at  all 
times. 

Let  the  maximum  value  of  grid  coordinate  1  be  R,  and  let  (n+2)At  be 
the  new  time  level  at  which  boundary  values  of  ?  or  0  must  be  given.  We 


write  IrQ+1  -  r0**  \ 

rH-2  „  n+1  _  IZK-LJ _ s  At 

*K-l,j  5K-1.J  Ax^  SK-l,jat 

(Ua+1  -  o0*1  ) 

n+2  n+1  *°K.j  .  A„ 

°K,J  •«kj  '  ^  Vl.i4' 


(i) 


(ii) 
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Both  schemes  (i)  and  (11)  were  tested  and  found  to  yield  comparable 
results.  The  discussion  will  be  confined  to  shceme  (11)  where  the 
boundary  value  of  the  velocity  Is  extrapolated  from  the  Interior  field 
value  at  the  earlier  time  step. 

Figure  12(a)  shows  the  wave  profile  at  28  seconds  that  resulted 
when  the  input  boundary  £  was  set  equal  to  zero  for  24  seconds.  A 

wave  pocket  has  been  caught  in  Its  parrage  toward  the  right  hand  boundary. 
In  Figure  12(b)  at  72  seconds,  the  wave  packet  has  vanished  through  the 
open  right  hand  boundary  simulated  by  means  of  equation  (11)  above. 


pig.  12  Wave  Profiles  for  Deep  Water  Configuration  (h  *  150  ft.) 

With  Open  Downstream  Boundary.  £  (1, j)  *  0,t£24  sec. 

Ax  ■  16  ft,  At  ■  0.025  sec,  n  *  0.026,  p  *  1 

6.  Wide  Area  Grid — Breakwater  Parallel  to  Beach 

As  a  test  of  the  ability  of  the  model  to  predict  longshore  effects, 

a  plane  wave  with  parameters  the  same  as  those  considered  above,  was  _ 

allowed  to  travel  past  a  breakwater  extending  half  way  across  the 

1904  x  1840  ft  (580.7  x  561.2  m)  grid  selected,  the  bottom  contour 

elsewhere  corresponding  to  the  beach  configuration  already  described 
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(h  *  ISO  to  -3  ft.).  Grid  dimensions  were  120  x  116,  with  Ax  ■  Ay  -  16 

ft  (4.9  m).  The  vertically  uniform  breakwater,  16  ft  thick,  was 

positioned  in  the  row  1  -  93,  where  the  water  depth  is  30.7  ft  (9.4  m), 

and  extended  928  ft  (283  n)  from  the  edge  of  the  grid,  parallel  to  y. 

It  was  thus  about  360  ft.  (110  m)  offshore.  Lateral  boundaries  to  the 

grid  were  closed.  A  constant  Manning  coefficient  of  0.026  was  assumed. 

There  was  no  artificial  smoothing. 

Figure  13  shows  schematically  the  positions  of  the  wave  crests  at 

t  *  64  seconds.  It  was  obtained  by  printing  symbolic  depths,  h  ,  such 

s 

that  £  >  0,  h  »  1;  0,  h  “0,  and  linking  the  0’s  by  shading, 

s  s 

Hence  white  areas  correspond  to  water  above  z  ■  0,  and  shaded  areas  to 
water  below  this  level.  Clearly,  there  is  a  diffraction  of  the  wave 
fronts  around  the  breakwater.  Substantial  values  of  V  occurred  in  the 
sheltered  region,  also  on  the  seaward  side  of  the  breakwater. 


Fig.  13  Wave  Crests  at  64  Seconds  for  the  Beach-Breakwater  Config¬ 
uration 

Ax  mAy  ■  16  ft  ,  At  ■  0.2  sec  ,  n  *  0.026 
7.  Conclusions 

A  model  solving  vertically  integrated  three-dimensional  equations 
with  borrowings  from  lineary  wave  theory,  has  been  shown  to  account  for 
short  waves  in  deep  water  very  well.  Problems  due  to  numerical  turbu- 
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lence  are  likely  to  arise,  however,  in  shallow  water  where  the  small 
amplitude  approximation  no  longer  holds,  and— -most  significantly — wave 
breaking  leads  to  a  multi-valued  solution  for  the  water  level  at  a 
point.  This  difficluty  can  be  overcome  if  the  shallow  zone  is  treated 
as  a  buffer  of  very  high  bottom  friction,  and  some  internal  energy  loss 
(smoothing)  is  employed.  In  effect,  a  beach  is  replaced  by  a  shock- 
absorbent  wall.  Such  an  approach  is  likely  to  be  satisfactory  for 
computations  where  offshore  processes  rather  than  beach  processes  are  of 
interest. 
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Appendix  1  -  The  Finite-Difference  Equations 
of  the  General  Wave  Model 

(See  Section  4  for  a  description  of  symbols) 


Cj  .i  ~  (Hu°  tj1^1  -  Hu11  u11"^  ^i+1/2 

1>j  1>J  +  (Hui+l,j  Ui+l,j  Hui,j  °if j  1+1/2 
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The  time-varying  land/water  boundary  for  the  beach  model  was 
handled  by  a  method  described  in  McHugh  (1976).  Immediately  following 
the  second,  fourth,  sixth,  etc.,  time  steps,  a  call  is  made  to  a  program 
subroutine  that  determines  whether  flooding  of  a  dry  square  can  occur, 
and  if  so,  transfers  water  into  the  dry  square.  The  criterion  for 
flooding  is,  in  the  case  of  a  one-dimensional  bottom  slope  parallel  to 

X:  St-1  > 

The  transferred  water  volume  Av  is  calculated  from  the  formula 
Av  =■  0.52b  (q_^j+  hi>:j)3/2  At 

where  b  *  3.33  is  the  Francis  weir  constant  (  Giles  (1962)).  This 
formula  is  obtained  by  time  integration  of  a  variable  crest  weir  flow 
over  the  interval  2 At,  assuming  the  water  level  in  the  flooding  square 
exceeds  the  bottom  level  in  the  flooded  square  during  all  of  2At.  The 
argument  for  this  assumption  in  the  case  of  the  beach  model  is  as 
follows:  First,  define  a  critical  depth  H£  below  which  a  newly  flooded 
square  is  considered  still  "dry".  This  strategy  was  adopted  to  avoid 
certain  numerical  problems  (McHugh  1976).  In  the  present  study  “  0.1 
ft.  (0.03  m).  With  wave  period  8  seconds  and  amplitude  1  ft.  the  water 
only  rises  approximately  0.01  ft.  in  0.02  seconds,  so  that  the  new  depth 
established  in  a  dry  square  after  flooding  during  time  2At  *  0.025 
seconds,  say,  will  be  very  much  less  than  Hft.  Consequently,  several 
applications  of  the  flooding  formula  will  be  necessary  to  establish  a  - 
new  depth  greater  than  Hc«  This  means  that  for  the  second  and  all 
subsequent  necessary  applications  c j  will  already  exceed  -h^  at 
the  beginning  of  the  interval  2At. 
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The  reverse  situation  to  flooding  of  a  grid  square  is  ebbing  of 
water  from  it.  In  general,  when  ebbing  occurs  a  negative  depth  will  be 
established  by  the  computer  program.  This  represents  a  volume  debt  to 
the  wet  field  that  must  be  restored  by  the  program  to  conserve  water 
volume  (  McHugh  1976). 


•  / 
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Fig.  6  Quasi-Beach  Wave  Profiles  at  72  seconds;  At  •  0.0125  sec 
0.026  (a)  Ax  -  16  ft.,  (b)  Ax  -  16,  8,  4  ft.  ’ 
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Fig.  8  Quasi-Beach  Wave  Profiles  at  72  and  73  Seconds  With  Coabined 
High  Bottom  Friction  and  Artificial  Snoothing. 

Ax  -  16,8,4,2  ft..  At  -  0.025  sec.,  n  -  0.026,  h  >  20  ft.;  n« 

0.2,  14  ft.  <  h  <  20  ft.,*  o  -  5,  h  <  14  ft. 
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Fig.  9  Beach  Wave  Profile  at  76  Seconds 
Ax  -  16  ft..  At  -  0.0125  sec.,  n  •  0.026 
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Fig.  10  Beach  Wave  Profiles  at  76  Seconds 

Ax  ■  16, 8, A,  At  *  0.025  sec.,  o  *  0.026,  h  >  20  ft.;  n  »  0.2,  10 
ft.  <  h  <  20  ft.;  n  0.5,  h  <  10  ft. 
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Fig.  13  Wave  Crests  at  64  Seconds  for  the  Beach-Breakwater  Config¬ 
uration 

Ax  ■  Ay  ■  16  ft  ,  At  «  0,2.  sec.,  n  ■  0.026 


IV.  PREDICTION  OF  PENETRATION 


It  is  well  known  that  the  physical  properties  of  bottom  sediments 
in  the  coastal  zone  are  highly  variable  both  in  space  and  time.  The 
nature  of  the  bottom  sediments;  mud,  silt,  sand,  or  cemented  sediments, 
can  have  a  large  effect  on  coastal  processes,  particularly  surface 
waves.  The  purpose  of  this  effort  was  to  support  the  Naval  Coastal 
Systems  Center  in  the  development  of  a  technique  for  rapidly  measuring 
marine  sediment  properties.  One  approach  that  has  potential  is  the 
instrumented  penetrometer  (I,  2,  3,  4).  The  objective  of  the  research 
was  to  review  existing  penetrometer  theory  and  instuments  and  indicate 
it's  potential  for  meeting  the  needs  for  classifying  bottom  sediment 
types. 

Recently  approaches  have  been  proposed  in  which  soil  strength 
properties  can  be  obtained  directly  from  continuous  measurement  of  cone 
and  sleeve  friction.  In  1973,  (5)  experiments  were  performed  in  which 
projectiles  were  penetrated  into  kaolin  clay  targets  with  known  and 
controlled  properties.  In  their  experiments  a  variety  of  projectile 
diameters  and  weights  as  well  as  nose  shapes  were  tested,  and  it  was 
found  that  for  kaolin  clay  the  side  wall  resistance  is  a  significant 
portion  of  total  resistance  and  that  at  low  velocities  the  shape  of  the 
nose  had  no  effect  on  the  magnitude  of  side  wall  resistance.  In  1973, 
(6)  a  dynamic  penetrometer  was  instrumented  with  decelerometer,  cone  and 
sleeve  load  cells.  Penetration  into  different  type  of  soils  was  inves¬ 
tigated.  This  investigation  indicated  that  such  an  instrument  is 


0 


capable  of  providing  shear  strength  characteristics  of  marine  sediments 
The  results  of  these  experiments  was  that  failure  mode  during  Impact 
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penetration  is  basically  of  a  static  nature  and  cavitation  does  not 
occur  for  impact  velocities  of  less  than  19  f/s,  and  that  the  friction 
ratio  concept  developed  from  static  cone  penetration  tests  for  assessing 
soil  types  can  be  extended  to  impact  penetration  tests. 

The  sediment  properties  which  have  been  related  to  penetration 
include  sediment  type  (1,  3),  shear  strength  (2,  4,  5,  6,  7),  sensi¬ 
tivity  (2,  4),  density  (2,  4,  5,  6),  viscosity  (5),  shear  rate  effects 
(6),  and  friction  angle  (7).  Typically  shear  strength  values  or  pro¬ 
files  are  determined  from  either  deceleration  (2,  3,  4)  or  tip  and 
sleeve  resistence  (6,  7). 

From  a  review  of  present  state-of-the-art  is  feasible  to  construct 
an  instrumented  penetrometer  to  measure  bottom  sediment  shear  strength 
and  classify  bottom  sediments.  The  instrumentation  would  consist  of  an 
accelerometer  with  a  range  of  _+25  g,  a  tip  resistence  gauge  with  a  range 
up  to  1000  Kg  and  a  frictional  sleeve  gauge  with  a  range  up  to  750  Kg. 
The  calibration  of  the  penetrometer  would  initially  be  done  on  labora¬ 
tory  samples  and  then  field  tested  on  a  variety  of  sediments.  The  cost 
for  a  research  instrument  would  be  in  the  range  of  $5,000  to  $6,000,  an 
operation  system  would  be  less  expensive,  but  a  lowest  cost  figure  can't 
be  determined  at  this  time. 
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PREDICTING  SURF  ZONE  CONDITIONS  FROM  AN  OFFSHORE  BUOY  SYSTEM 
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ABSTRACT 

A  buoy  system  and  two  surf  prediction  models  are 
being  developed,  tested,  and  evaluated  for  obtain¬ 
ing  wave  breaker  statistics  and  longshore  currents. 
Input  data  Is  based  upon  measurements  from  an  off¬ 
shore  location  outside  the  surf  zone. 

Data  from  pressure-transducer  wave  and  tide  gage, 
a  two-axes  electromagnetic  current  meter,  and 
temperature  sensors  are  digitized  and  sent  In  a 
serial-digital  format  to  an  off-line  processor. 

Both  models  used  to  describe  the  nearshore  waves 
and  longshore  current  velocity  are  based  upon 
linear  theory.  Results  of  measurements  Indicate 
that  prediction  of  breaker  heights  within  an 
accuracy  of  10  to  IS  percent  Is  feasible  using 
offshore  measurements. 

INTRODUCTION 

The  coastal  environment  plays  a  key  role  In  many 
of  the  Navy's  vital  operations.  Of  specific 
Interest  to  the  Navy  are  such  parameters  as  near¬ 
shore  waves,  tides,  currents,  and  surf  conditions. 
Actual  measurements  of  these  parameters  within  the 
active  surf  zone  Is  a  difficult  task,  particularly 
If  It  Is  to  be  done  In  a  tactical  situation.  The 
surf  zone  Is  a  very  hostile  environment  on  most 
occasions.  And,  It  Is  where  the  surf  Is  hostile 
that  Information  concerning  its  statistics  are 
most  needed.  Measurements  outside  the  surf  zone 
have  been  done  more  extensively  and  with  con¬ 
siderably  more  success.  Also,  data  processing 
techniques  have  been  developed  for  offshore 
conditions  as  have  environmental  Instrumentation. 
Therefore,  the  use  of  these  developed  processes 
and  hardware  coupled  together  with  an  effective 
reliable  prediction,  can  provide  valuable  infor¬ 
mation  to  the  coastal  investigator. 


This  paper  will  report  on  the  efforts  of  the  Naval 
Coastal  Systems  Center  (NCSC)  to  develop  an  off¬ 
shore  buoy  system  for  measuring  wave  and  current 
parameters  and  related  prediction  models  that  take 
these  parameters  and  generate  surf  zone  statistics, 
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The  specific  details  of  the  buoy  system.  Its  appli¬ 
cation,  the  software  algorithms,  the  data  processl,.? , 
and  test  results  are  given  by  Pldgeon. 1 

DATA  BUOY  SYSTEM 

The  design  of  the  data  buoy  system  Involved  con¬ 
siderations  of  many  factors.  One  of  the  major 
objectives  of  the  design  was  to  develop  a  measure¬ 
ment  system  that  would  operate  In  shallow  water 
outside  of  the  surf  region  and  would  be  able  to 
measure  the  hydrodynamic  parameters  needed  as  off¬ 
shore  Inputs  into  the  prediction  model  for  deter¬ 
mining  the  statistics  of  wave  breakers  and  the 
longshore  currents  Induced  by  these  breakers.  A 
secondary  consideration  was  the  development  of  a 
system  that  might  either  lead  to  or  provide  signif¬ 
icant  information  for  the  design  of  an  operational 
buoy  system  that  could  be  used  in  Fleet  operations 
and  provide  outputs  detailing  the  near-surf  and 
surf  regions  of  an  area  under  consideration.  Al¬ 
though  it  waa  realized  that  this  would  be  a  multi¬ 
phase  program,  ultimate  end  products  were  consid¬ 
ered  along  the  various  development  stages  of  the 
present  system.  While  considerable  program  devel¬ 
opment  was  done  using  the  Hewlett  Packard  HP-9825 
programmable  calculator,  additional  processing  was 
done  using  microcomputers  in  preparation  for  adap¬ 
ting  the  buoy  system  to  perform  on-board  data 
processing  utilizing  state-of-the-art  micropro¬ 
cessing  components  within  the  buoy. 

The  buoy  system  measures  waves,  tides,  two  or¬ 
thogonal  components  (horizontal)  of  current,  water 
temperature.  Internal  temperature,  and  system 
calibration.  Waves  and  tides  are  simultaneously 
measured  by  a  Statham,  Model  UC-3,  pressure  cell. 

A  Marsh-McBlrney,  Model  711,  tvo-axls  electro¬ 
magnetic  Induction  current  meter  measures  water 
motion  In  the  horizontal  plane.  Water  temperature 
is  measured  both  Inside  and  outside  the  unit  by 
using  two  Analog  Devices,  Model  AD590,  two-terminal 
Integrated  circuit  temperature  transducers.  In 
addition,  a  leak  detector  circuit  was  added  to  the 
buoy  to  detect  the  presence  of  water  In  the  bottom 
of  the  unit.  All  signal  outputs  from  the  system 
Bensors  are  analog  voltages  end  are  tied  Into  a 
16-channel  analog  to  digital  multiplexer  having 
12-blt  accuracy  prior  to  transmission  to  shore. 

The  sensors  and  all  associated  electronics  are 
contained  in  a  cylindrical  pressure  housing  capable 
of  withstanding  pressures  equivalent  to  30  meters  of 
water  depth.  Weighing  less  than  25  kilograms  In 
air,  the  buoy  Is  easily  deployed  by  divers  and  Is 
normally  attached  to  a  stand  that  has  been  previ¬ 
ously  jetted  Into  the  bottom  for  support. 
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A  major  consideration  In  the  design  of  the  buoy 
system  was  the  method  of  getting  the  analog  signals 
from  the  sensors  to  the  processing  equipment. 
Desirably,  the  data  would  be  sent  In  digital  form 
in  order  to  obtain  the  most  immunity  from  noise 
and  other  degrading  factors.  Also,  In  order  to 
avoid  having  many  signal  lines  between  the  buoy 
and  the  receiving  station,  a  serial  not  a  parallel 
system  was  needed.  A  very  important  consideration 
In  the  design  of  this  Interfacing  system  was  the 
formatting  and  makeup  of  the  data.  With  future 
on-board  processing  In  mind,  a  format  that  would 
be  directly  compatible  with  mocrocomputers  and 
other  similar  processors  is  what  is  needed.  Each 
of  these  considerations  fit  very  well  Into  an 
interfacing  system  that  has  been  developed  by  Dynex 
for  another  NCSC  task.  This  system,  called  Infor¬ 
mation  Transmission  System  (ITS)2,  takes  multiple 
analog  data  channels  (up  to  128),  multiplexes  them 
into  a  12-blt  digitizer  (giving  one  millivolt 
resolution  out  of  four  volts),  converts  the  par¬ 
allel  output  data  Into  two  elght-blt  bytes,  which 
are  sequentially  applied  to  a  standard  UART  (Uni¬ 
versal  Asynchronous  Receiver  Transmitter) ,  from 
which  the  data  are  sent  in  a  bit-serial,  digital 
format.  The  serial  format  Is  completely  compatible 
with  computer  standards  for  serial  data  transfer, 
l.e.,  meets  RS-232C  standards.  Therefore,  any  com¬ 
puter  processor  with  a  standard  RS-232  serial  in¬ 
terface  can  accept  the  data  from  this  buoy  system 
without  any  further  modification. 

An  Important  consideration  In  the  transfer  of  data 
from  the  measuring  source  to  the  processing  system 
is  the  preservation  of  the  data  as  nearly  as  mea¬ 
sured  as  practical.  Secondly,  the  ability  to  Input 
directly  Into  the  processor  with  minimal  reformat¬ 
ting  means  more  efficiency,  fewer  errors,  and  less 
processing  costs.  For  on-board  processing,  both  of 
these  considerations  are  almost  mandatory.  These 
factors  are  the  prime  ones  considered  In  the 
development  of  the  ITS  system,  and  they  also  made 
It  a  prime  candidate  for  the  data  buoy  system.  The 
addition  of  ITS  to  the  buoy  system  allowed  consid¬ 
erable  flexibility  in  how  the  data  could  be  handled, 
once  It  was  obtained  In  the  buoy.  The  data  could 
be  transmitted  via  standard  telemetry  links  with 
relatively  low  bandwldths.  Or  the  data  could  be 
directly  connected  to  a  beach  site  via  cable.  This 
latter  version  was  selected  for  these  preliminary 
tests  alnce  both  data  signals  and  buoy  power  could 
be  cabled  between  the  beach  site  and  the  buoy  with 
a  modest  cable  having  as  few  as  four  conductors — 
two  for  power  and  two  for  signals. 


Since  It  was  desirable  to  record  data  at  selected 
Intervals,  it  was  decided  to  Input  the  data  signals 
directly  from  the  buoy  cable  into  a  standard  Bell 
103  telephone  modem,  having  an  auto-answer  capa¬ 
bility.  The  RS-232  data  format  at  a  300  baud  rate 
Is  completely  compatible  with  this  type  modem. 

Thus,  when  It  was  desired  to  take  data  or  to  cheek 
the  buoy  operation,  all  that  one  had  to  do  was  to 
call  up  the  buoy  from  HCSC  or  the  Dynex  office  with 
either  an  acoustic  coupler  (which  was  used  by  NCSC) 
or  a  modem  (which  was  used  by  Dynex) .  When  the 


modem  at  the  beach  tower  auto-answered,  the  hand¬ 
shake  was  completed  at  the  originating  end.  At 
this  point,  data  began  to  stream  in  and  continued 
until  the  originating  source  turned  off  Its  lock¬ 
up  tone  or  the  calling  phone  was  hung-up.  The 
data  that  was  received  could  then  be  used  In 
several  ways.  It  could  be  checked  to  ascertain 
proper  buoy  operation,  it  could  be  used  to  deter¬ 
mine  the  real-time  conditions  of  the  Gulf  of  Mexico 
at  that  time,  or  It  could  be  directly  recorded  or 
input  Into  the  processor  for  data  processing. 

In  order  to  be  able  to  sample  the  wave  data  rapidly 
enough  to  detail  with  accuracy  the  wave  profile, 
the  wave  data  were  sampled  most  often  of  the  data. 
This  was  done  by  putting  the  wave  data  on  every 
other  channel.  Interdispersed  beween  the  wave  data 
were  the  other  signals.  With  the  ITS  system  run¬ 
ning  at  300  baud  wave  data  were  sampled  almost 
five  times  per  second.  This  was  sufficiently  fast 
enough  to  adequately  determine  wave  shape,  (It  Is 
necessary  In  the  data  processing  to  determine  the 
peak  and  trough  of  each  Individual  wave  as  well  as 
Its  period.  Therefore,  an  accurate  reproduction 
of  the  wave  profile  Is  necessary.) 


PREDICTIVE  MODELS 

Using  the  wave  buoy  to  provide  a  time  history  of 
hydrodynamic  pressure  and  the  two  horizontal  com¬ 
ponents  of  water  motion  at  a  point  seaward  of  the 
surf  zone,  a  scheme  was  developed  to  compute  the 
following  parameters: 

a)  Average  breaker  height  (H)  and  maxlmun 
breaker  height  (H  ) 

IB 

b)  Average  breaker  position  (X)  and  maximum 

breaker  position  (X  ) 
m 

c)  Breaker  type  (l.e.,  plunging,  spilling  or 
surging) 

d)  Average  longshore  current  (Vj)  and  maximum 
longshore  current  (V1bj) 

Wave  theory  used  for  this  study  Is  linear  wave 
theory,  This  theory  Is  able  to  treat  waves  in  in¬ 
termediate  and  shallow  water  and  is  easily  program- 
able  onto  a  microcomputer.  It  la,  of  course, 
limited  to  waves  of  small  height  to  wave  length 
ratios  (low  steepness),  and  small  wave  height  to 
water  depth  ratios.  In  addition  to  wave  and 
current  parameters,  bottom  slope  (P) , elevation  of 
buoy  off  the  bottom  (dl)  and  orientation  of  current 
meter  (y-axla,  parallel  to  shoreline  and  x-axla, 
normal  to  shoreline)  are  also  required  Input  para¬ 
meters  (see  Figure  1).  Data  la  sampled  at  a 
constant  time  Interval,  At,  over  an  Interval  of 
time,  t,  to  produce  N  data  samples,  where  N  equals 
t/ht.  Usually  t  will  be  In  the  range  of  five  to 
20  minutes  and  At  will  be  0.25  to  one  second. 


60 


PRELIMINARY  COMPUTATIONS 

The  preliminary  confutations  Include  determining 
the  mean  preaaure  (P),  the  preaeure  change  between 
wave  creet  and  wave  through  for  each  wave  (AP;  J-1 
•  •  • M) ,  the  wave  period  for  each  wave  (T) ,  the  Bean 
velocity  coaponente  (V  and  U)  and  the  wave  crest 
current  coaponente  (U^  and  ) . 

The  wave  paraaeters  coaputed  for  the  location  of 
the  lnstruaent  are: 


water  depth 
and  for  each  wave: 
wave  height 

wave  length 

Wave  angle 


—  +  dl 
Pg 


„  .  AP.1  COSH  (2ffh/Lp 

J  Pg  COSH  (2irdi/L^) 

i1  -  TANH  (2i»h/L.) 

TJ  Lj  3 

8  -  TAN-1 

9J  TAN  (uCJ) 


SIMPLIFIED  MODEL-PREDICTION  EQUATIONS 

For  the  prediction  of  breaker  heights,  the  simpli¬ 
fied  model  utilizes  pressure  measurements  (uncor¬ 
rected  for  the  depth)  from  the  bottom  mounted 
pressure  transducer  and  computes  a  breaker  height 
from  the  following  equation: 

>i  «r],/i 


«b 


■  [• 
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Pg 


+  h. 


This  computation  Is  made  for  each  Individually 
measured  wave  recorded  at  the  data  buoy. 

The  predictive  model  also  gives  a  prediction  of  the 
maximum  longshore  current  from  the  measured  average 
longshore  current  at  the  buoy  and  the  predicted 
breaker  depth  obtained  from  the  mean  breaker  height 
as  given  by: 


(y 


1 )  Breaker  Depth 


h^  -  H^  (1-4.78  +  128  }»  8  <  .15 


2)  Breaker  Angle 

sin  sin  e 


hj  C 


J 


J 


CJ  ■  ^  •  cbj  ■ 

3)  Breaker  Wavelength 

“  CbJ  *  T1 

4)  Breaker  Height 

VHJ  KrJ  *  K.J  *  Kf  j 
where 


Kf J  ‘  -98  •  KrJ  ‘ 
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The  set  of  breaker  values  are  computed  for  j  * 
1 • • *M.  For  this  the  maximum  breaker  height 
can  be  Identified,  furthermore 


u  _  1 


r  _  1 


M 


?b  "  «  A  8»J 

The  positions  of  h^  and  hus  can  be  located  on  the 
nearshore  chart. 


DETAILED  MODEL-PREDICTION  EQUATIONS 

The  prediction  of  the  surf  zone  parameters  Involves 
a  set  of  equations  which  must  be  applied  simulta¬ 
neously.  The  equations  use  the  lnstruaent  site 
wave  paraaeters  ,  T ^ ,  L^ ,  and  8^  to  predict  the 

same  paraaeters  at  the  break  point  H^,  Tb>  L^,  8^. 
It  la  assuaad  that  the  period  of  each  wavs  reaalns 


constant  so  that  T,  *  T, 
follows  (3):  3 


bj’ 


The  equations  are  as 


PREDICTION  PROCEDURE-DETAILED  MODEL 

The  simultaneous  solution  of  the  set  of  equations 
needed  to  make  a  prediction  la  too  complex  for  a 
practical  scheae.  Rather,  the  solution  can  be 
approximated  by  an  Iteration  scheae  In  which 
Initial  values  for  selected  paraaeters  are  computed 
from  an  Input  data.  The  formula  for  the  Initial 
value  of  the  wave  height  Is 

H*  -  .95  h1/5  H  4/5 
DJ  J 


for  tine-delayed  neaaureaents  to  the  breaker  region 
for  a  similar  group  of  breakers. 


With  this  Initial  value,  the  set  of  equations  In 
the  surf  prediction  model  can  be  computed  explic¬ 
it  el v  in  order,  that  Is,  finding  h^,  use  this  to 

fl»_  and  finally  H^.  The  difference  between 

and  H*bj  should  be  small  (  a  few  percent)  and, 

If  not,  R.  .  can  be  used  In  place  of  H*  In  a 
second  Iteration.  J 


BREAKER  TYRE 

The  prediction  of  breaker  type  Is  based  upon 
breaker  wave  characteristics  and  the  beach  slope. 
The  empirical  formula  suggested  Is: 

*^>1  >1.2  SPILLING  BREAKER 

V8 

<1.2  PLUNGING  BREAKER 


SURGING  BREAKER 


LONGSHORE  CURRENT  MODEL 

The  prediction  of  longshore  currents  Is  based  upon 
the  theory  of  linear  waves  on  a  plane  sloping  beach, 
The  formula  for  the  maximum  longshore  current 

and  the  longshore  current  at  the  average  point 
are  given  by  (4): 

Vm  -  20.7-B(gHb)5i  SIN  20b 


and 


V 
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Two  different  approaches  to  the  verification  of 
the  predicted  helghte  were  utilised.  The  first 
approach  developed  was  a  technique  using  a  transit 
with  vertical  graduations.  This  vertical  scale 
was  calibrated  for  known  locations  (anchored  floats) 
within  the  breaking  sone.  Vertical  displacement  of 
the  float  was  measured,  giving  a  reasonably  accurate 
measurement  of  waves  immediately  prior  to  breaking. 

A  comparison  of  measured  versus  predicted  breaker 
heights  for  the  predictive  models  Is  presented  In 
Table  1.  Although  It  was  not  possible  to  make  a 
one-for-one  comparison  between  the  predicted  breaker 
heights  and  the  measured  breaker  heights,  their 
respective  statistics  provided  adequate  data  for 
comparison. 

It  can  be  observed  that  the  detailed  model  over¬ 
predicts  the  heights  in  all  cases,  whereas  the 
simplified  model  comes  closer  to  the  measured  con¬ 
ditions,  In  some  cases  under-predicting,  but  In 
others  over-predicting.  However,  on  the  average, 
the  simplified  model  Is  very  close  except  for  the 
largest  of  the  heights.  For  this  test  date  the 
average  wave  periods  are  running  close  to  six 
seconds.  It  should  be  noted  that  periods  of  less 
than  three  seconds  are  not  Included  due  to  fre¬ 
quency  cut-off  of  pressure  transducer. 

A  brief  examination  of  the  effect  of  beach  slope 
on  the  results  given  for  the  detailed  model  was 
done.  A  slope  of  1:20  was  used  for  the  data  pro¬ 
cessing.  At  this  point,  this  value  Is  felt  to  be 
too  large.  1:30  to  1:50  probably  fits  the  con¬ 
ditions  more  closely.  Using  the  more  gradual  slope 
In  the  detailed  model,  the  prediction  results  for 
the  cases  tested  was  less,  but  not  by  as  much  as 
the  results  for  this  test  date  Indicate.  The  slope, 
as  well  as  other  factors,  need  to  be  further 
Investigated  for  model  refinement. 


FIELD  TESTS  AND  RESULTS 

To  date  the  buoy  has  been  deployed  on  five  differ¬ 
ent  occasions  to  acquire  data  for  various  sea  state 
conditions.  Prior  to  each  deployment  and  upon 
return  to  NCSC,  the  entire  system  undergoes  a  check 
procedure  to  verify  operability  and  calibration. 
Problems  encountered  with  the  use  of  a  cable  be¬ 
tween  buoy  site  and  shore  station  have  been  the 
most  prevalent.  Strong  longshore  currents  associ¬ 
ated  with  the  passage  of  storm  fronts  have  caused 
the  cable  to  break  or  to  pull  away  from  the  buoy 
on  several  different  occasions. 


To  properly  evaluate  the  predictive  models,  It  was 
desired  to  compare  actual  measurements  of  breaker 
heights  on  an  Individual  basis  with  corresponding 
predicted  heights.  This  necessitates  time  correl¬ 
ation  of  recorded  buoy  measurements  with  observed/ 
measured  breaker  measurements.  Due  to  the  disper¬ 
sive  nature  of  waves  approaching  the  breaking  sone, 
the  sequence  of  waves  arriving  at  the  buoy  will  not 
necessarily  be  the  same  sequence  that  arrive  at  the 
breaker  sone.  However,  In  general,  the  measure¬ 
ments  of  a  group  of  waves  can  be  compared  with  that 


It  should  be  noted  at  this  point,  that  the  results 
of  the  comparisons  of  the  models  with  observations 
are  very  encouraging.  The  results  indicate  that 
for  breaker  heights  up  to  two  meters,  the  detailed 
model  predicts  breaker  heights  to  an  accuracy  of 
the  order  of  ten  percent.  For  the  higher  energy 
waves,  the  simplified  model  does  almost  as  well, 
but  does  less  well  for  the  shorter  wave  periods 
with  lower  amplitudes. 

The  second  approach  to  model  verification  Involves 
the  use  of  wave  data  collected  from  the  U,  S,  Army 
Coastal  Engineering  Research  Center  Field  Research 
Facility  at  Duck,  North  Carolina.  Breaking  wave 
data  recorded  from  a  series  of  wave  gages  positioned 
throughout  the  active  surf  tone  will  be  compared 
with  predicted  breaker  heights  whose  input  data 
have  been  simultaneously  recorded  from  gages 
positioned  just  outside  the  surf  tone.  These 
detailed  time  correlated  measurements  will  provide 
valuable  calibration  data  for  both  predictive  models. 
Results  of  this  experiment  will  be  published  at  a 
later  date. 
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CONCLUSIONS 

Regarding  the  over-all  buoy  system  and  lta  opera¬ 
tion,  the  results  of  thla  study  have  demonstrated 
that  thla  system  la  a  very  reliable  one,  has  the 
ability  to  measure  very  accurately  the  required 
wave,  tide,  and  current  Inputs  for  data  processing, 
Is  easily  deployable  by  divers  from  a  small  boat 
with  minimum  equipment  required,  and  Is  readily 
adaptable  to  a  variety  of  Fleet  operations  and 
applications.  It  has  demonstrated  that  It  can 
withstand  the  hostile  environment  of  the  near-surf 
zone  within  the  limits  tested  and  has  the  ability 
to  transmit  lta  data  via  a  wide  variety  of  methods 
to  the  end  user,  thus  making  this  buoy  system  a 
prime  candidate  for  Fleet  operations. 

The  models  presented  in  this  study  provide  a  means 
for  making  predictions  of  the  desired  surf  para¬ 
meters.  Each  model  Is  limited  In  application 
because  of  the  assumptions  contained  in  the  under¬ 
lying  wave  theory.  The  models  do  not  account  for 
wave  Induced  pressures  near  the  break  point,  result¬ 
ing  from  wave  crest  accelerations.  Furthermore, 
the  bottom  topography  Is  assumed  to  be  planar  and 
therefore  real  beach  profiles  and  longshore  bars 
cannot  be  treated.  Finally,  the  longshore  current 
predictions  neglect  temporal  variations  in  current 
speed  and  spatial  variations  associated  with  near¬ 
shore  topography. 


TABLE  1.  COMPARISON  OF  MEASURED  VERSUS 
PREDICTS  BREAKER  HEIGHTS  FOR 
SIMPLIfILD  AND  DETAILED  MODELS 
TEST  DATE  7  DECEMBER  1979 


TIME 

AVERAGE 

BREAKER 

HEIGHT 

. 

«H 

II 

«B  V 

4J  T3 

V  O 

• 

(0 

w  x 

O  X 

X 

11:15-11:18 

.69 

.83 

11:18-11:22 

.70 

.81 

11:22-11:26 

.76 

.85 

AVERAGE 

.72 

.82 

.73 

11:15-11:26 

11:28-11:31 

.69 

.83 

11:31-11:35 

.79 

.86 

11:35-11:39 

.75 

.84 

AVERAGE 

11:28-11:39 

.74 

.84 

.76 

AVERAGE  ALL 

TIMES 

11:15-11:39 

.73 

.83 

.75 

It  is  concluded  that  the  software  developed  for  the 
support  operations  of  the  buoy  system  is  efficient 
in  terms  of  processing  capabilities  of  speed,  size, 
and  complexity,  can  be  adapted  to  on-board  buoy 
processing  without  considerable  modification,  Is 
flexible  and  readily  adaptable  to  a  wide  range  of 
user  applications,  and  the  study  results  further 
show  the  software  to  be  reliable  and  accurate. 

It  is  further  concluded  that  both  of  the  prediction 
models  used  In  this  study  give  very  good  results, 
often  with  prediction  errors  of  only  a  few  percent, 
that  the  detailed  model  allows  for  a  greater  flexi¬ 
bility  and  range  of  environmental  inputs  but  does 
require  more  complexity  In  processing,  and  that  the 
simplified  model,  although  limited,  gave  accuracies 
that  are  well  within  the  requirements-  of  many  user 
applications. 
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FIGURE  1.  Schematic  Diagram  of  Beach 
and  Definition  of  Variables 
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